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Abstract 



The Zeldovich approximation is applied to study the accretion of hot and cold 
dark matter onto moving long strings. It is assumed that such defects carry 
a substantial amount of small-scale structure, thereby acting gravitationally 
as a Newtonian line source whose effects dominate the velocity perturbations. 
Analytical expressions for the turn-around surfaces are derived and the mass 
inside of these surfaces is calculated. Estimates are given for the redshift 
dependence of the fraction of mass in nonlinear objects. Depending on 
parameters, it is possible to obtain = 1 at the present time. Even with 
hot dark matter, the first nonlinear filamentary structures form at a redshift 
close to 100, and there is sufficient nonlinear mass to explain the observed 
abundance of high redshift quasars and damped Lyman alpha systems. These 
results imply that moving strings with small-scale structure are the most 
efficient seeds to produce massive nonlinear objects in the cosmic string model. 
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I. INTRODUCTION 



The cosmic string model may be the most promising alternative to inflation-inspired 
theories of structure formation (see Ref. for recent reviews). Many predictions of the 
model, for example a scale-invariant primordial spectrum of density perturbations and 
the amplitude of cosmic microwave background (CMB) anisotropics on large angular scales 
0], are at least qualitatively in agreement with observations. 

The main problem facing the cosmic string model is the computational difficulty in 
analyzing its predictions. Cosmic strings will form during a phase transition in the very 
early Universe provided that certain well known topological criteria @] are satisfied by the 
matter theory. It is also well established that the network of strings, once formed, will 
approach a "scaling" solution i.e. the distribution of strings will become independent of 
time if all lengths are scaled to the Hubble radius. 

However, the exact form of the scaling solution is not known. At any time t, there are 
two components to the string network: loops with radius R <^t, and "long" string segments 
with curvature radius Rc ~ t. The best available numerical simulations [^] indicate that 
more mass is in the long strings than in loops. A further uncertainty pertains to the amount 
of small-scale structure which builds up on the long string segments as a consequence of 
the involved dynamics of the string network. Long straight strings with no small-scale 
structure carry no local gravitational potential. Their only effect is to lead to a conical 
structure of space perpendicular to the string 0]. Their transverse velocities will typically 
be close to the speed of light since the underlying dynamics is governed by the relativistic 
wave equation. Strings carrying small-scale structure, on the other hand, generate a local 
Newtonian potential in addition to the conical distortion of space. 

Given the lack of knowledge about the details of the initial string distribution, it is 
important to study all mechanisms by which strings can give rise to structures. Loops 
will produce roughly spherical objects by gravitational accretion p|, long strings without 
small-scale structure generate planar overdensities [^|10[) while long strings with small-scale 
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structure induce - since they act as Newtonian gravitational line sources - filamentary objects 



IMll- 



In this paper we will study filament formation, i.e. the accretion of both cold dark matter 
(CDM) and hot dark matter (HDM) onto a moving Newtonian line source, by means of the 
Zeldovich approximation and its suitable adaptation to HDM |TB[. The Zeldovich ap- 



proximation is a Lagrangian perturbation technique which is an improvement over standard 
linear perturbation theory and is widely applied in cosmology. 

The accretion of CDM onto static string loops was analyzed a long time ago |]l6l . Subse- 



quently, it was realized [|T3;|lMl ^^^^ - contrast to theories based on adiabatic perturbations 
- the cosmic string model also appears viable if the dark matter is hot. The effects of string 
loops on HDM were then analyzed carefully by solving the collisionless Boltzmann 



equation, i.e. keeping track of the entire phase space distribution of the dark matter par- 
ticles. Studies of nonlinear clustering about static string loops were also performed pO[ by 



means of N-body simulations. At the time when string loops are produced, they typically 
have relativistic velocities. Hence, it is important to study the accretion of dark matter by 
moving seeds. Bertschinger applied the Zeldovich approximation to study the case of 
CDM. He found that the total nonlinear mass is to a first approximation independent of the 
velocity Ui of the loop. The decrease in the width of the nonlinear region is almost exactly 
compensated by the increase in the length of the structure. The study of clustering onto 
moving string loops was recently p2[ generalized to HDM. In this case, the total nonlinear 



mass was shown to decrease with increasing Ui. Nevertheless, it was demonstrated that 
sufficient accretion to form nonlinear objects by a redshift of ^ = 4 takes place. 

Since the string network simulations indicate that long string would be more important 
for structure formation than loops, attention turned to the formation of string wakes. The 
accretion of CDM was analyzed in detail by means of the Zeldovich approximation in Ref. 
p3| . In Ref. [|I^ , a modification of the Zeldovich approximation was introduced which made 
it possible to study the clustering of HDM in a simple way. This approximation was tested 
against the more rigorous approach obtained by solving for the phase space distribution of 



dark matter particles by means of the collisionless Boltzmann equation. The two methods 
were shown to give very similar results. 

As emphasized by Carter it is likely that long cosmic strings have a substantial 
amount of small-scale structure, and that hence the Newtonian line potential gives rise to 
more important gravitational effects on dark matter particles than the velocity perturbation 
induced by the conical structure of space perpendicular to the string. The accretion of 
cold [|l3l and hot dark matter onto a moving string with a wiggle has been analyzed 



several years ago. However, a better description of the coarse-grained small-scale structure 



is by an effective Newtonian gravitational line source P,|lT||12[. Hence, in this paper we 
will study the accretion of CDM and HDM onto a moving Newtonian line source by means 
of the Zeldovich and modified Zeldovich approximations, respectively. We will neglect the 
velocity perturbation which is due to the conical structure of space. The main result is 
that even in the case of HDM, nonlinear structures form at a redshift close to 100. For 



static line sources, our results reduce to those of Ref. [^]. As we shall see, for reasonable 
values of the parameters of the cosmic string model, either with CDM or HDM, the theory 
is compatible with the present observational constraints coming from the abundance of high 



redshift quasars. Some of our results have been obtained previously in |T1| and [O. However, 



our analysis is more detailed and goes further, and some of the conclusions concerning the 
accretion of hot dark matter are different. 

The outline of this paper is as follows. In the following section we give a brief review of 
the standard and modified Zeldovich approximations. In Section 3 we study the accretion 
of cold dark matter onto a moving line source, and in Section 4 we modify the analysis for 
hot dark matter. Section 5 concerns an estimate of the total nonlinear mass as a function of 
redshift z for the cosmic string filament model. We conclude with a discussion of the main 
results. We work in the context of an expanding Friedman- Robert son- Walker Universe with 
scale factor a{t). Newton's constant is denoted by G, h denotes the Hubble expansion 
parameter in units of 100 km s^^ Mpc"^, and Q is the ratio of energy density to critical 
density. 



II. THE ZELDOVICH APPROXIMATION 



The Zeldovich approximation is a first order Lagrangian perturbation method based 



on analyzing the trajectory of individual particles in the presence of an external gravitational 
source. In an expanding Universe, the physical position of a dark matter particle is given 
by 

r(q,t) = a(t)[q-V'(q,t)] , (1) 

where q is the comoving coordinate and if) is the comoving displacement due to the gravi- 
tational perturbation. 

By combining the Newtonian gravitational force equation, the Poisson equation for the 
Newtonian gravitational potential, using conservation of matter, and expanding to first order 
in '0 one arrives at the following equation (valid in a matter dominated Universe) for the 



comoving displacement of a background particle p^,pi|,p2 

d'^ih 2adib 3a , ,, 

where a dot means total time derivative and S stands for the source term. 

Up to this point, the analysis is quite general. We now specialize to the case in which 
the perturbations are generated by a long straight cosmic string with small scale structure 
whose strength is given by the value of GX, where X = fi — T, fi and T being the energy 
per unit length and the tension of the string, respectively. If the position of the string is 
denoted by r'{t) = a(t)q'(t), then the source term 5* is given by: 

^ 2XG q-q' 



|q — q'p 



(3) 



In general, there are two sources of perturbations induced by strings, the wakes produced 
by the conical structure of space perpendicular to the string, and the Newtonian gravitational 
line source given by (|^). The effect of wakes can be modeled by a nonvanishing initial 
velocity perturbations towards the plane behind the string. However, since we are interested 
in exploring filament formation, we shall neglect wakes. Hence, the appropriate initial 



conditions for the perturbed displacement are limits) = dxl){ts)ldt = 0, where ts is the time 
when the perturbation is set up (see below). Using the Green function method, the solution 
of Equation (|^) may be written in the following form 

V^(q,t)=A(t)I(q,t)-D2(t)J(q,t) , (4) 

with 

Here, Di and D2 are the two independent solutions of the homogeneous equation, and W 
is the Wronskian. We take Di and D2 to be the growing and decaying modes, respectively. 
In the case of a spatially flat (Einstein - de Sitter) Universe (to which we will restrict our 
attention) and for ti > t^q we have a = (t/to)^''^, Di = a and D2 = a~^^^. 

For cold dark matter tg is the true initial time (the time when the string sets up 
the perturbation), while for hot dark matter we must take neutrino free streaming into 
account. In principle, one should study the evolution of the phase space density by means 
of the Boltzmann equation and integrate over momenta to obtain the density distribution. 



However, as demonstrated in |TB] in the case of planar accretion, it is a good approximation to 
use Equation (H) but with modified initial conditions. Since the effective perturbations due to 
the seed survive only on scales greater than the free-streaming length [T^ Aj ~ v{t){z{t) + l)t, 
where v{t) is the r.m.s. velocity of the hot dark matter particles at time t, we can take for 
HDM a scale-dependent starting time tg, namely the time when |q — q'| = Xj{ts{q)) (where 
|q — q'l is the relative distance between the seed and the dark particle) which yields: 

where Aq = Viti\/ Zi + 1. The implementation of this initial condition in the basic solu- 
tions, in order to account for the effects of free streaming, is called the modified Zeldovich 



approximation |15|. In the case of static filaments, the above HDM starting condition is 



unambiguous. However, for moving filaments tliere is an ambiguity. For a fixed q, the value 
of ts{q) changes as the filament moves. In fact, ts(q) may initially be smaller than t and 
later, when the filament is closer to q, increase to a value larger than t. We will specify how 
to deal with this ambiguity in Section 4. 

A comment is in order concerning the motion of the line source. If we suppose the string 
source is straight (and infinite), pointing along the z axis, then the acceleration along the 
z axis does not affect the motion of the string. That is, we can interpret q' in Equation 
(0) as a vector position in the perpendicular plane to the string. The same is also valid for 
the motion of a fluid dark particle, since its 2;-motion is not perturbed by the string source. 
Thus, in the formulae related to the filament all vectors are to be thought as belonging to 
the xy plane. For instance, in Equation (|^) we have to take q = {qx, Qy) ■ 

The dark matter particles are initially moving away from the string with the Hubble flow. 
The gravitational action of the string slows this expansion. At the time when dx/dt = the 
particles decouple from the Hubble flow. This is called "turn-around". At any time t, we 
can compute the value of which is turning around. We denote this quantity by qni,x- The 
above turn-around condition, together with Equation ([l|) yields 

qnl,x = ^i'xiqnl, t) + tpx{(lnl, t) , (8) 
(X 

which determines the comoving turn-around surface. 
Using Equations (^), (|^) and (H) we get 

qni,x = 2alx{qni,t) + ^a~^/V^.(q„i, t) , 

= 2^x{<ini, t) + -a-^/'^Jx{qni, t). (9) 

In order to determine the turn-around surface we must integrate Equations and 
for which we have to specify the position of the string, i.e. c({t). We choose initial condition 
for the seed motion in such a way that q'(tj) = 0, and r'(tj) = Uj. Since the peculiar velocity 
of the string decreases as a(t)^\ we have 

(10) 

where qi = 3uiti/ai is the final comoving position of the string. 



q'(^) = qi 



1 - 



1/2 
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III. ACCRETION OF COLD DARK MATTER 



In this section we study the accretion of CDM onto a moving filament. By considering 
the motion of the string along the y direction (then q' = q\a)ey) and using ([T0|), 
and we obtain 



9tg \Gq^ 

5 ql + {.qy-qiY 



tan-^ ( ^y-^'^^A - tan 



ql + {qy - q, 

-I ( qy-qs \ 



0+2 

Ja: = ^AGg,.F(q,a), 
5 



(11) 
(12) 

where = qi{l — ^Jai/ag) denotes the position of the string at the time when clustering for 
HDM begins (for CDM Qs = 0), and Y stands for the integral 

a^da 



F(q,a) 



(13) 



qx + ky - q'ia)? 

We see that Y scales like a^/^ when a — >■ 00. In this case it is sufficient to keep only the 
leading term in a for the decaying mode: 



a-3/2j,. ^ ^XGa- 



(14) 



25 -^ql + {qy-q,y ' 

It is worth mentioning that the above equation follows from (0) and ([T3|) by neglecting 
terms of order of a^^^ (and lower). Within such an approximation we may neglect the a 
dependence of q' in the In and tan~^ terms in equation (|ll]), because its contribution is of 
the same order as terms neglected in (14). 

Using Equations (H)-(|l^ it is easy to shown that the turn-around condition (which can 
be evaluated at any redshift z) is given by 



1 

+ 10+ 



+2 



qy - qi 
qx 



tan 



-1 / qy - qi 
\ qx 



tan 



-1 ( qy-qs 
, qx 



(15) 



with 
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b(t) = ^XGa. (16) 



Note that the turn-around surface has the form of a tube whose perimeter is defined by (|T5[) 
and whose length is proportional to the comoving distance corresponding to the Hubble 
radius at time tj. Thus, the comoving length lc{ti) of the string at ti is given by 

UU) = ^{1 + ziU))U , (17) 

where ^ is a constant which determines the curvature radius of the long string network 
relative to t and whose value will be discussed in Section 6. For CDM the initial conditions 
are = 0, = with the right hand side of (|I5|) reducing asymptotically to 



tan- ^ -tan-M^ 



(18) 



where we defined Q'com ll + {% ~ li^"^ 

The perimeters of the turn-around surfaces at various redshifts are depicted in Figure 
1. Unlike in the static case, the string motion along the y axis gives rise to a suppression 
of the turn-around distance along the x axis. For moving loops, the same effect has been 



established by Bertschinger [^. It is easy to understand the shape of these curves. The 
filament starts its motion at |q'| =0 with a high velocity Mj. Due to the expansion of the 
Universe, the peculiar velocity decreases at later times as a{t)~^. As a consequence, the 
seed spends much less time in the region, say, qy < 0.9qi than in the region qy ~ gj. The 
integrated action of gravity is hence larger in the second region. This explains why the 
turn-around distance is greater near the final position of the filament (see also Figure 2). It 
should be noted, however, that if the filament is replaced by a string loop, the turn-around 
surface will be thinner nearby the end of the motion. This can be understood by the fact 
that the gravitational pull for loops {F ~ 1/r^) is less efficient at large distances than for 
filaments (F ~ 1/r). 

Figure 1 
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FIG. 1. The comoving perimeter of the turn-around surface for the moving filament, with initial 
velocity Ui = .2 and CDM. The curves are for the redshifts z = 200 (innermost contour), 50, 20, 5, 
3, and 1 (outermost contour). For z < 2.1 the turn-around consists of two contours. 

It is worth mentioning that the accretion is also suppressed (on small scales, qx ~ 0) near 
q = (0, qi), as happens with the small scale suppression due to free streaming. This is shown 
in Figure 2 where the turn-around redshift is plotted as a function of q^ for qy = qi/2, and for 
qy = qi where such a suppression is observed (dotted line and continuous line). These curves 
look like the case of a static filament in HDM models |25], however, here the suppression at 
a fixed time t ^ ti occurs only on scales much smaller than the suppression caused by the 
free streaming. This is readily seen by rewriting equation ([T8| ) for qy = qi 



qI = list + Kt) 



(19) 



where qx,st is the turn-around comoving position (the x component) in the case of a stationary 
filament. For small values of Ig^,!, the In term on the right hand side of (|l^) contributes with 
a very large negative value (the tan~^ terms are finite even for qy ^ qi), thereby implying 
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that turn-around is possible in that region just for very large times (compared to ti). This 
suppression is also a consequence of the peculiar motion of the string whose gravitational 
force on a cold dark particle acts in different (though correlated) directions at different times, 
which is not the case on large scales. Such a small scale suppression does not occur in the 
case of loops This can also be explained by taking into account the gravitational action 
of each defect. 

Figure 2 
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FIG. 2. Redshift versus Qx for different values of qy and Uj. 



The overall effect of the string motion is to distort the cross-section of the turn-around 
surface (a circle in the stationary case). However, the suppression along the x axis is more 
than compensated by the growth along the y axis (see Figure 1). In fact, the suppression 
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along X is only logarithmic in the string velocity (it goes as ln{l/qi)) while the length of the 
wake along the y axis increases linearly (proportional to gj). In contrast, the total amount 



of CDM accreted by a string loop is independent of the velocity (see Refs. p6| and 
Figure 3 shows the dependence of the maximum value qni,x{%) as a function of the velocity 
of the seed. 

Figure 3 
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FIG. 3. Variation of the maximum of qni,x{Qy) as a function of the velocity of the seed. The 
plots are for z = 1000 (lowermost curve), 200, 50, 20, 10, 5, 2.5, and 1 (uppermost curve). The 
difference for qni,xiQy, ^i) between Uj = and Ui = .2 falls from 25% for z = 1000 to 11% for z = 1. 

We have calculated numerically the total nonlinear mass (i. e. the mass inside the turn- 
around surface) as a function of the redshift. A reasonable approximation can be obtained 
by assuming that the area inside the turn-around curve defined by ([l^) is given by 



A — TlQcDM + Qi-^ y 



(20) 



where is the the comoving width of the turn-around volume at its thickest point (near 
% = qi). The second term in ( pO]) corresponds to the area of an isosceles triangle whose base 
is Ix and height is qi. This equation in fact yields a lower limit for the total mass that has 
gone nonlinear. Since a lower bound to Qni will be derived in Section 5, it is appropriate at 
this point to obtain a lower limit to the accreted mass. 
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The second term in (|20|) is proportional to Ui, the initial velocity of the string, whereas 
the first term in the limit Ui = reduces to the results of Ref. for the static case. 
Comparing the two terms, we see that the second term in (|20|), stemming from the motion 
of the string, gives the dominant contribution if 

(AG)f (^(t) + l)-i/W/2<l, (21) 
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where {XG)q is the value of AG in units of 10~^. For example, we see that if ti = teq, 
{XG)q = 0.5, h = 1/2, Q = 1 and Ui = 0.2, then the two terms are comparable at redshift 0, 
but the contribution originating from the second term in (|20|) , the contribution specifically 
due to the motion of the string, dominates for higher redshifts. 
The total mass accreted can be written as (see Equation (|17D) 

M{z, u^) = ipoU{zi + l)A{z, Ui) , (22) 

where {zi + 1) = {t^/tif'/'^ and A{z,Ui) is the area given by (pO|). As a consistency check, 
we note that taking the static limit (g^ = 0) in the above equations, the results of Ref. ||2^ 



are readily recovered. For the values of the parameters used above. Equation ( P^ yields: 
M ~ 5^ X W^MQA{z,Ui)lMpc^. 

To proceed further, it proves convenient to summarize the main properties of the turn- 
around curve described by (|l^). For CDM we have checked numerically: 

(i) The comoving area grows linearly with the scale factor (as a consequence of the fact 
that ln{a/ai) is a slowly varying function of a) for a ^ Oj. 

(ii) The area inside the turn-around curve grows almost linearly with the velocity of 
the filament Ui for Ui > .05 and is practically constant for Ui < .05 at redshifts under 
consideration (because in the latter case the term in (|20| ) independent of Ui dominates). 

(iii) Due to the motion of the filament, the accretion on large scales is strongly suppressed 
near the comoving initial position of the string (because the gravitational line source is far 
away for most of the time), whereas growth on small scales is suppressed at the end, when 
q' ~ Qi- 
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(iv) Using Ui = 0.2 and tj = teg we obtain at the redshift z = 4 a turn-around curve whose 
length is / ~ gj ~ 8h~'^Mpc and whose width in the thickest part is about 2h^^{XG)l^'^ Mpc. 
The area A inside the curve is approximately 16h~"^ {XG)l^^ Mpc^ which corresponds to a mass 
per unit length of M = ^Apot^gZeg - 50^ x 10^^ h~\\GYJ'^MQ. 



IV. ACCRETION OF HOT DARK MATTER 

Now we turn to the accretion of hot dark matter. The basic equations (|TT| - [l^) which 
describe accretion onto a moving Newtonian line source are applicable. 

We will deal with the ambiguity in the definition of ts{q) mentioned in Section 2 in the 
following way: Since we are only interested in the motion of dark matter in the direction 
perpendicular to the seed velocity, we may replace the initial condition by 

\'ix I 

This, however, seems to overestimate the free streaming effect, since it completely erases 
the accretion in the region near the initial location of the string (an improved analysis using 
the Boltzmann equation is in progress and will be reported elsewhere). In the following we 
will adopt ( P3[ ) as a generalization of the condition used for static seeds which is a good 
approximation for the real free streaming effect [|1^ . The quantities we are going to obtain 
in this case, such as the total mass accreted, are lower bounds to the exact ones. 

Since we must have ts > ti, condition ( ^3]) can be applied only if < Viti{zi + 1), 
otherwise the condition tg = ti is to be used even for HDM. 

Thus, only a few small changes must be made: In the turn-around conditions we substi- 
tute Qg = j^-^. Then for HDM, equation (^) may be cast as 



where = qi{l - \qx\/K), and Aj = Xo/^/a~= Xj(ti). 
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(24) 



As shown in Figure 4, the turn-around curve p^ ) is composed of two identical discon- 
nected closed curves, one for negative values of and the other for positive values of q^. It 
is apparent from the figure that the nonlinearities on small scales are highly suppressed by 
the free streaming condition - see (H). 



Figure 4 
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FIG. 4. The same as Figure 1, but now for HDM with velocity Vgq = .1. Contours (two sheets 
each one) are drawn for z = 20 (inner contour), 5, 1, and (outer contour). 



In order to compare to other cases let us rewrite (|2^) for qy = qi 

In 



ql + 



(25) 



,<ll + <li\<ix\\i VJ ' 

Recalling that equation ( pSf ) applies for \qx\ < Xi (for larger values of the CDM result 
is valid) we see that the first In term on the right hand side of (^) is negative while the 
second one is positive. At early times, when ql « qf those two terms reduce to ln{\qx\/Xi) 
and a numerical analysis shows that for Vi < Ui the ratio of this term to q^ cdm very 
small. This means that for a rapidly moving filament the suppression of the turn-around 
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near Qy — qi is governed by the seed motion and is the same for both CDM and HDM. 
For Vi > Ui the suppression due to the free streaming is important also for qy = q^. These 
features can be seen from Figures 2 and 5, graphs which plot the redshift corresponding 
to the onset of nonhnearity as a function of qx for different velocities of the string. For 
instance, for the value Ui — 0.2, the curves for qy — q^ in the cases of CDM (Figure 2) and 
HDM (Figure 5) are essentially the same, whereas for the value Ui = 0.05, the HDM curve 
is substantially lower at small values of q^- Figure 5 also shows that there are no nonlinear 
structures for redshifts greater than about 90. 

Figure 5 
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FIG. 5. Same as figure 2, now for HDM with velocity Vgq = .1. 
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We know that free streaming completely erases the nonlinearities on scales smaller than 
the free streaming length. Equation ( p^ shows this behavior as a consequence of the term 
ln{ql/Xj), which is independent of Qy and for small \qx\ is highly negative implying that 
turn-around can only occur at late times. As we have seen in the previous section, for 
Qy < Qi the tan~^ terms are also important in determining the turn-around on large scales 
(they are finite even for — > 0). Then, the presence of the new tan~^ term in (|25| ) will 
not significantly change the results, and it follows that for a moving filament with HDM, 
the wake is strongly erased in the thinnest region near the starting position of the seed (see 
Figure 4). 

The total nonlinear mass accreted in the HDM model can be estimated by considering 
all the mass inside and between the two disjoint curves (see Figure 4). We then approximate 
the turn-around curve by a trapezoid. Such a figure figure is obtained when the height of 
the isosceles triangle considered in the CDM case is reduced to /y, thereby accounting for the 
free streaming suppression. By comparing Figures 1 and 4 it can be seen that the difference 
in the nonlinear mass between HDM and CDM is large at high redshift and decreases at 
later times. The length ly of the turn-around curve along the y axis depends on the redshift 
and tends to for very late times (e.g. for 2; < 2 if tj = teq)- According to Figure 4, a simple 
(underestimating) approximation is to take ly = qi — qy{K), where qy{\i) is given by the 
intersection between the line q^ = —\qy/qi + \ and the approximated turn-around curve 
for CDM (see the comments below (|20D ). The result is ly = qi — Xiqi/{Xi + Ixf^)- Hence, the 
area of the nonlinear region can be estimated as 

where QIjj^m stands for the right hand side of (pSj). The first term is analogous to the one 
introduced in the CDM case, the first term in equation (^), which represents the accretion 
onto a static filament (at the comoving position qi) with HDM. For t sufficiently greater 
than ti the second term in ( |2^ ) dominates if the same condition for CDM, Equation (PT]), is 
satisfied. In such a case we can parameterize the result as follows: 
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MHDMit,ti,Ui) = f{t,ti,Ui)McDM{t,ti,Ui), (27) 

where f{t,ti,Ui) is the function which parameterizes the mass suppression in the case of 
HDM. For large t, the result with HDM asymptotically approaches the CDM result. Hence 

limt^oofit,ti,Ui) = l. (28) 

The numerical analysis shows that for the cosmological and cosmic string parameters used 
(see Section 3), at a redshift of z = 10, f(t,teq,Ui) ~ 0.5, and for 2 = 4, f{t,teq,Ui) ~ 0.8. 
The function f(t,ti,Ui) = 1 — (^x +i' /2)'^ used in (^) furnishes the following approximate 
values: / = 0.5, / = 0.85, and / = 0.96 respectively for z = 10, z = 4 and z = 0. We have 
taken = 2|ga;| with given by (^) (see following discussion). 

V. THE NONLINEAR DENSITY PARAMETER Qnl 

To illustrate the results of the previous sections, let us now compute the fraction Qni of 
the critical density accreted onto filaments both for CDM and HDM models. Naturally, by 
considering only the accretion onto filamentary structures we will be underestimating the 
total contribution to Qni due to the cosmic string network. Let us first consider the case 
where the initial string velocities Ui and times t are such that the second term in (|20|) , the 
term specific to the motion of the string, dominates. This is the case of primary interest to 
us, since in the other limit (g^ = 0) the results for static filaments are recovered (see eqs.(20) 
and (26)). 

As argued earlier (see Section 3), to estimate flni in the case of moving filaments with 
CDM, one may approximate the turn-around curve by an isosceles triangle whose height 
and base are respectively ly = qi and l^, where is the width of the thickest region of the 
wake. A good approximation for the width is 2\qx\, where is given by ([T9|) without the 
term ln[ql/{ql + qf)]. This term leads to a suppression of small scales, whereas for large 
scales it is small in comparison to ln{a/ai). In this way, by considering the mass of the wake 
for < qy < qi, such a term may be safely neglected. The resulting expression for q^ is then 
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the same as for a static filament located at qi since the initial time tj. Also, the condition 
that the string is not near qi for all times after tj, can be effectively described replacing ti 
by the rescaled time te// = ct^j? where < a < 1. The time tg// can be thought as a sort 
of initial time for which the Newtonian line starts to act effectively as a static source for 
accretion near qi. The value of a may be fixed by comparing the exact value of q^ obtained 
numerically with the approximation given by qx,static{t,<yti). As can be easily seen, for the 
values of the parameters considered here one has a ~ 10^^. In the limit a ^ 1, the relative 
error for sufficiently large times may be written as ln{a)/ln{t/ti), and for ti = teq, it is about 
15% at the present time (see Figure 3 for comparison). 

The total mass accreted onto the filament during an arbitrary time interval At = t — ti 
{ti < t < to), is the mass contained in the Lagrangian volume inside the turn-around surface. 
Since the length of the turn- around tube is proportional to the Hubble radius at time ti, the 
accreted mass is given by (see Section 3) 

-<-^'^€/iKs)"'faR' (-) 

This is the total mass that has gone nonlinear at time t > ti onto a Newtonian line formed 
at time tj. Note that the nonlinear mass grows as t^^^ and not as t^/^ as in the linear theory 
for static seeds. The difference comes from the fact that the length along the direction 
of string motion depends only on the velocity of the seed, and is, to the first order of 
approximation, time independent. As expected, for slowly moving strings, the first term in 
(pop dominates, and the usual linear perturbation theory growth of the total nonlinear mass 
is recovered. This can also be seen by using ( pT)) to determine the lower cutoff velocity when 
(p9|) is applicable, and substituting in (^9]) for the velocity Ui. This yields the extra factor 

of {Z{t) + 1)^1/2. 

From the scaling solution of the string network we know that the number of long strings 
per physical volume at time t created in the interval between ti and ti + dti assume the form 

n{t,ti)dti = ^dti, (30) 
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where z/ is a constant which gives the number of long string segments passing through a 
volume at any given time (this number is independent of time as a consequence of string 
scaling). Using this number density we obtain the fraction of the critical density of objects 
formed by accretion of CDM onto moving filaments in the interval between t^g and t as 

Qni{t) = QT^Gt^ f n{t, U)M{t, U)dti , (31) 

where M{t,ti) is given by (^). 

Then by direct integration of (|3TD we obtain to a first approximation 

1/3 



'ISAG 



— -M- - 



In 

2^i>iu,)o.2{XG)l/'h{z{t) + ir'/\ (32) 



where ('Uj)o.2 is the value of Ui in units of 0.2. It follows that for Ui = 0.2 and {XG)q = 0.5 
the value of Qni{to) at the present time to is about 2^z/. Note that the above result is valid 
for values of Ui for which (pID is satisfied. For smaller initial string velocities, by solving 
(pT|) for Ui and inserting in (0)), we obtain in the static limit 

^niAt) = y^AGz> (^pj ' In ^ ^z>|^(AG)6(z(t) + ly^h^ . (33) 

For HDM we have approximated the turn-around curve by a trapezoid (see the discussion 
below equation (^)). The velocity depending area Ah dm inside such a curve is given by 

AhDM = fit, t,, U,)AcDM = ( 1 " Jy^-^fjif) ' ^^'^^ 

As mentioned before, this area includes the holes inside and between the curves of Figure 
4. This means that all In terms in (|25| ) are neglected. As a matter of fact, the suppression 
terms due to free streaming are important only for small scales which have been by the above 



approximation included in Ah]:,m- Besides, the term in (p5D arising from the motion of the 
string may also be neglected by the same reason as discussed earlier. As a consequence, the 
expression of for HDM reduces to the same one derived in the CDM model. 
The value of in the case of HDM can be estimated by following (|3TD : 
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^ni""^' = evrCt' / n{t, ti)f{t, ti, Ui)McDM{t, Ui)dti . (35) 
Since the integral is dominated at the lower end, we can approximate the result as 

where, as we have seen above, f{t, t^g, Ui) grows from / ~ at t = teg to / ~ 0.96 at t = tg- 
For the sake of completeness, let us now write the contributions from the first terms in 
equations (pOD and (|2^) to the fraction of nonlinear mass. As discussed earlier, these terms 
correspond to the fraction of the accreted mass onto static filaments. Following the previous 
procedure we get 

^Sf - ^4<-^G { ln{±-)^ (37) 



5 yO^^eq / y^^eq ^ 



(compare with (^)) and for HDM 



riLst 



where 1 — x is the suppression factor due to free streaming, x = xit^teq) falls from x = 1 
for t = teg to X ~ at t = to- Notice that Qni,st varies with {z + l)^-*^ while the contribution 
due to the string motion, equations (|3^ ) and (|36|), goes as {z + 1)~^/^. This fact explains 



why Qni,st is important just for small z, when compared to the terms depending on the seed 
velocity. 

These results show that even with hot dark matter (given the string parameters used 
above), accretion onto filaments can easily lead to Qni = 1 at the present time. 

VI. CONCLUSIONS 

We have studied the accretion of hot and cold dark matter onto moving gravitational line 
sources and applied the results to the cosmic string filament model of structure formation. 
Our main result is that even for initial relativistic motion of the strings, nonlinear structures 

20 



can develop at redshifts close to 100 if the dark matter is hot (a point not realized in earlier 
work |11|). Thus, string filaments are much more efficient at clustering hot dark matter at 
early times than string wakes. 

The total fraction Qni of nonlinear mass depends quite sensitively on the string parame- 
ters used, i.e. on the initial string velocity Ui, on the strength AG of the Newtonian potential, 
and on the values of the constants ^ and 9 which determine the curvature radius and number 
of strings in the scaling solution. However, the uncertainties partially cancel. For example, 
the product uC, is a measure of the total mass density in strings and can in principle be 
determined from numerical simulations of cosmic string evolution. Also, as AG increases, 
the value of Ui will decrease, but in this case the cancelation is only partial since Ui cannot 
exceed 1 and hence in the limit AG the filamentary accretion will disappear (and the 
only effect of the strings will be the wakes caused by the velocity perturbations). 

According to the numerical simulations of Allen and Shellard @] (which, however, are 
based on strings obeying the Nambu-Goto action and hence not directly applicable to strings 
with a large amount of small-scale structure), the value of is about 15 in the radiation 
dominated epoch and slightly less during matter domination. Taking the value (AG)6 = 0.5 
which corresponds to the largest possible effect of small-scale structure for GUT scale strings, 
we find that even for hot dark matter, all of the mass in the Universe has gone nonlinear by 
a redshift of about 5. This implies that the cosmic string model - provided the small-scale 
structure dominates - is easily compatible with the cosmological constraints from both high 



redshift quasar and damped Lyman alpha system p8[ abundances. 

It is interesting to compare the results derived here with those from the cosmic string 
wake [1^ and loop p2[ scenarios. If the dark matter is hot, wakes do not produce any 



nonlinearities until a redshift of close to 1. Loops do produce nonlinear objects at high 
redshifts, but objects of mass IO^^Mq only form at a redshift of about 4. Hence, we conclude 
that string filaments are the most efficient way of forming high redshift massive nonlinear 
structures in the cosmic string model. Unfortunately, it is not known how much small- 
scale structure cosmic strings are endowed with. The cosmic string evolution simulations 
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at present do not have the required range to resolve this issue, and analytical analyses are 
needed to address this question. In light of the sensitive dependence of high redshift structure 
formation to the amount of small-scale structure, it is crucial to improve our understanding 
of the small-scale structure on strings. 
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